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In disordered superconductors, the local pairing field fluctuates in space, leading to the smearing 
of the BCS peak in the density of states and appearance of the subgap tail states. We analyze the 
universal mesoscopic contributions to these effects and show that they are enhanced by the Coulomb 
repulsion. In the vicinity of the quantum critical point, where superconductivity is suppressed by 
the "fermionic mechanism" , strong smearing of the peak due to mesoscopic fluctuations is predicted. 
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Superconductive (s-wave) state is characterized by a 
gap A in the quasiparticle spectrum and the coherence 
peak (CP) in the density of states (DOS) above the gap, 
p{E) = on ReE/\/E 2 — A 2 . According to classical re- 
sults 0, IH , impurity scattering does not affect this pic- 
ture, as long as the time-reversal invariancc (TRI) is 
not broken. Yet, a number of experiments demonstrate 
considerable suppression of the CP and appearance of 
subgap (E < A) states with the increase of disorder 0- 
jHJ ■ A mechanism leading to gap smearing without invok- 
ing any TRI breaking was proposed 40 years ago in the 
seminal paper @ (see also Ref.0). It was shown that the 
effect of a (phenomenologically introduced) short-scale 
disorder in the Cooper attraction constant, A = A+<5A(r), 
is formally equivalent to the one produced by magnetic 
impurities Another mechanism Q relates smearing 
of the CP with a finite inelastic lifetime of quasiparti- 
cles this effect becomes exponentially weak at low 
temperatures, T <C T c . Finally, it has been recently 
demonstrated that an apparent DOS smearing seen in 
tunnelling experiments may be due to electric fluctua- 
tions in the environment 

In recent few years an upsurge of interest in experimen- 
tal studies of strongly disordered (non-granular) super- 
conductors has been seen, evidenced, e. g., by Refs. \12\ - 
16j . Two basic classes of these materials distinguished 
by the value of the electron concentration are known 
(for a review, see Ref. [17|). Below we will focus on 



strongly disordered superconductors with high (typical 
metallic) electron density and strong Coulomb interac- 
tion fl3L UBI . 16 1 , where the fermionic mechanism of su- 
perconductivity suppression by disorder [18( is operating. 

In this Letter we show that mesoscopic conductance 
fluctuations [l9[ provide a universal lower bound for the 
DOS smearing effects in any disordered superconduct- 
ing thin films, effective down to T — 0. For the case 
of thin films (thickness d is below the low-temperature 
coherence length £o)j the strength of this smearing is 
completely controlled by the film dimensionless conduc- 
tance g = 2irh/e 2 Ro ~> 1, and the critical conductance, 
g c = ln 2 (?i/T c oT«)/27r, for the fermionic mechanism of su- 



perconductivity suppression [l8(. Here T c q is the transi- 
tion temperature in the clean system (we put Boltzmann 
constant fc# = 1) and r* = max{r, r(d/Z) 2 }, where r is 
the elastic scattering time and I = vpr is the mean free 
path. 

We find that the average DOS schematically shown in 
Fig.Q]is characterized by two energy scales: The width T 
measures the broadening of the BCS peak, while r ta ;i de- 
termines the exponential decay rate of the subgap DOS, 
(p{E)) oc cxp{-[(i; g - £)/T tail ] 3 / 2 }. The tail in the av- 
eraged DOS is a manifestation of the local gap inhomo- 
geneity due to randomness in impurities' configuration. 
In the zero-temperature limit, r ta ii and T r ta ii are 
given by 
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The most important feature of the result ((T) is a sharp 
growth of the DOS broadening in the vicinity of the 
quantum critical point, g = g c . The same dimension- 
less parameter, 5d ~ 1/3(3 — 9c), is known to control the 
disorder-induced smearing of the thermal transition [20| . 

Our quantitative analysis is presented below. 

Mean-field structure of the superconducting state. — We 
start the analysis of the superconducting state at the 
mean-field level, working with disorder-averaged quanti- 
ties and neglecting spatial fluctuations of the order pa- 
rameter. The effect of the Coulomb interaction on the 
properties of disordered superconducting films is usually 
treated in terms of an energy-dependent Cooper ampli- 
tude A(C), with £ = ln(l/£Y*) being the logarithm of the 
running energy scale (hereafter we set h = 1). In the case 
of the screened Coulomb interaction, A(£) obeys the RG 
equation fl8j : 



dX/dC = X 2 - Xl Xl = l/2ng, 



(2) 



where we neglected the triplet sector contribution and 
conductance renormalization assuming £/ ' g 1 . Equa- 
tion ((2) describes a competition between Cooper instabil- 
ity and Coulomb suppression of Cooper attraction, with 
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FIG. 1: (color online). Schematic view of the average DOS in 
a dirty superconducting film (solid line). Broadening of the 
BCS peak (dashed line) is mainly described by the semiclassi- 
cal approximation (dotted line), with the full DOS containing 
a significant tail of the subgap states. 



the initial condition A(0) = Ao at E ~ l . The RG flow 
© drives A(C) to infinity at C* = (2A s )ln[(A + A s )/(A - 
X g )], which signals the superconducting transition with 
the critical temperature T c = T^ 1 e^^* [l8|: 




(3) 



Another important feature introduced by the Coulomb 
interaction is the energy dependence of the pairing po- 
tential A(e) [2^. The latter is defined through the qua- 
siclassical Gorkov function in the Matsubara representa- 
tion, F(e) = A(e)/[e 2 +A 2 (e)] 1 / 2 . The function A (e) can 
be obtained from the self-consistency equation (SCE) 



A(e) 



(4) 



where e is the fermionic Matsubara energy, and the 
energy-dependent Cooper amplitude is given by 



A(e, e 1 ) = Ao - A 2 ln[l/ max(e, e')t* 



(5) 



A description based on Eqs. ((4]) and ((5j), where all ener- 
gies are retained, is alternative to successive elimination 
of high-energy degrees of freedom by the RG evolution of 
A(C). The logarithmic correction to Ao in Eq. (O corre- 
sponds to the last term in RG equation $2$. At the same 
time, the first term in Eq. ([2]) is automatically taken into 
account by the summation over energy in Eq. (j4|). An 
approach based on Eqs. ((!]) and ([5]) is equivalent to the 
RG in the determination of T c , but can be used also to 
find A(e) in the superconducting phase at T < T c . 

Equation (j4|) generalizes the SCE in the BCS theory. 
To find the energy dependence of A(e) at large energies, 
e T c , we rewrite Eq. 0) as a linear integral equation 
in terms of the logarithmic variable £, which is readily 



solved by reducing to a differential equation owing to a 
simple form of the kernel ([5]). As a result, we arrive at 



A(e) = A(T)[(e/T c )^ + (e/T c )-^]/2, 



(0) 



valid for £ > T c . Here A(T) is a function of tempera- 
ture which should be determined from the full equation 
(HJ), where the region of small energies, e ~ T c , becomes 
important. Neglecting a slow e-dependence of A(e), we 
conclude that A(T) is related to T c as in the standard 
BCS theory [in particular, A(0) = 1.76TJ. 

According to Eq. ([6]), high-energy (E 3> T c ) electrons 
experience a larger value of the effective pairing potential 
[201 ] . This effect is most pronounced in the limit of strong 
suppression of superconductivity, T c <C T c q, when the 
overall enhancement becomes large: A(t j ,T 1 ) ^> A(T). 

Smearing by inhomogeneities. — The mean-field theory 
developed above describes disorder-averaged quantities. 
In the presence of a quenched disorder, the order pa- 
rameter becomes non-uniform and the sharp BCS peak 
gets broadened. Analytical description of this effect is 
complicated due to the failure of the perturbation theory 
at E — > A. Earlier experience @, [2l[ suggests that the 
problem can be conveniently tackled in two steps: 

1. First one has to find the correlation function 
(A(r)A(r')). Since this is a thermodynamic quan- 
tity involving contributions from many energies, it 
can be obtained in the Matsubara representation 
by a regular perturbation theory. 

2. Then behavior of electrons in the field of a spatially 
fluctuating A(r) can be considered independently 
for each real energy E. 

This scheme based on the perturbation theory is appli- 
cable provided that disorder smearing is small. 

In dirty superconductors, diffusive motion of electrons 
is described by the Usadel equation [221 ] supplemented by 
the SCE (j4l^ Various types of disorder, such as magnetic 
impurities[8|, [23| or fluctuating coupling constant 0, [2l[ , 
can be easily incorporated into the scheme as random 
fields in the Usadel equation. The situation with uni- 
versal mesoscopic disorder (intrinsic fluctuations of the 
potential disorder) we are considering is different: Since 
the Usadel equation is already written for the ensemble- 
averaged quantities, mesoscopic potential disorder can- 
not be included there as some extra fluctuating field. To 
find the correlation function (A(r)A(r')), one has to go 
beyond the Usadel equation, consider two replicas of the 
system and average over soft diffusive modes [20I ] , similar 
to calculation of the universal conductance fluctuations 
(UCF) 0. 

Mesoscopic fluctuations of the order parameter. — In 
the presence of disorder, the SCE (j4|) contains two sources 
of disorder: mesoscopic fluctuations of the coupling con- 
stant, Adi s (e,e';r) (in its Coulomb part), and meso- 
scopic fluctuations of the Gorkov function, Fdi s (e;r). 
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These quantities exhibit fluctuations even for a uni- 
form order parameter since they are governed by dif- 
fusive motion of electrons sensitive to mesoscopic dis- 
order. Due to the SCE, this fluctuations will result in 
an inhomogeneous contribution to the order parameter: 
A(e; r) = A(e) + Ai(e; r). The latter, in turn, will mod- 
ify F which therefore can be represented in the form 
F = F + (<9F /<9A)Ai + F dis . Substituting this in the 
SCE (jl} and linearizing we get an equation for Ai(e) in 
the Fourier representation: 

A l( e, q )-.T£A( £ , e 0g^Ai( £ ',q) 

<9A(e',q) 

= ttTJ2 [X(e,e')F dis (e',q) + A dis (e, e', q)F (e')] ■ (7) 
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FIG. 2: (color online). Plot of the function K' 1 (T/T c ), where 
K(t) = L 2 (t)N{t) is defined by Eqs. © and JTTJ. 



In the 2D case, the main contribution to the CP smearing 
comes from large scales [see Eq. (fTB"])] . much exceeding the 
correlation length r c ~ £ ~ ' D/T c for mesoscopic fluc- 
tuations of Fdis and Ad; s - Therefore it suffices to consider 
Ai(e) at zero momentum which will be implied below. 

Inverting the matrix in the left-hand side of Eq. (0, 
we express Ai(e) in terms of Fdi s and Adis- To study 
the CP smearing we need the small-energy limit of Ai(e) 
with e ~ T c , where the effect of fluctuations is to modify 
A(T) -> A(T) + 6A(T), with 6A(T) given by 



A A (T\ T ( T \ ( 27TT )' 

5A(T) = L (-)— ) 



E1,E 2 >0 

x [A(£i,e 2 )Fdi s (£2) + Adis(ei,£2)F (e 2 )]. (8) 

Here Lq is the fluctuation propagator at zero momentum 
and frequency in the BCS theory: 



V(|)=^E 



A 2 (T) 

23 ' 



<S = v/e 2 + A 2 (T), (9) 



with the asymptotic behavior L ps 47r 2 F 2 /7C(3)A 2 (T) 
at T ^ T c , and L = 1 at T = 0. 

Therefore, smearing of the CP and behavior near the 
gap edge are determined by a single number, /(0) = 
(SASA)q—o, which can be easily obtained from the corre- 
lation functions of Fdi s and Adis with the help of Eq. (jHJ . 
Since the Coulomb correction to A already contains a 
closed loop (return probability) , the correlation functions 
(FdisFiis), (F dis A dis ) and (A d isA d is) are given by the one-, 
two- and three-loop diagrams in soft diffusive modes, re- 
spectively. The overall contribution is given by [HI 



irDA(T) / T 
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where K(T/T C ) = L 2 {T /T C )N(T /T c ) (see Fig.©, and 

A(T) 

ei,e 2 >0 



N 
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11) 



with N ps 14C(3)A(T)/7r 3 F at T — > T c , and TV = 1 at 
T = 0. Replacement of l/.g 2 by l/g(g - g c ) in Eq. (JTUJ 
is due to high-energy contributions, T c < e < r^T 1 , where 
A(e) is enhanced according to Eq. ([6]). In the limit of 
strong T c suppression, g — g c < g Cl the leading source of 
disorder comes from mesoscopic fluctuations of the return 
probability in (AdisAdis)- Equation (|10p is consistent with 
our previous result in the vicinity of T c [20j , generalizing 
it to arbitrary temperatures T <T C . 

Mean-field density of states. — The average DOS, 
(p{E)) = po Re(cos(9(F, r)), is expressed in terms of 
the spectral angle 8 which satisfies the Usadel equation, 
(D/2)W 2 9 + iEsm9 + A(r) cos6> = 0, with a random or- 
der parameter, A(r) = A + SA(r) [here A = A(T)]. 
Integrating out short-range degrees of freedom one get s 
an equation for the long-range behavior of 9{E, r) (3.l2l|: 



D 



V 2 9 + iE sin9 + A(r) cos 9 ~ A ncos9sm9 = 0, (12) 



where the depairing strength is expressed in terms of the 
correlation function /(q) = (5ASA)q as 



'/ 



2 f f(q) d 2 q 



A J Dq 2 (2tt) 5 



(13) 



This expression has been originally derived in Ref. [6| 
for the 3D geometry. The last term in Eq. (fl~2"|) coin- 
cides with the depairing term due to magnetic impurities 
derived by Abrikosov and Gorkov (AG) Q. It leads to 
the broadening of the CP (shown by the dotted line in 
Fig. [I]), with the hard gap at E^ G = (1 - ?; 2/3 ) 3/2 A . 

In the marginal 2D case, mapping to the problem of 
magnetic impurities should be done with care. Con- 
trary to the 3D geometry, now the integral in Eq. (fT3")) 
is logarithmically divergent at small q. An appropri- 
ate cutoff can be established by retaining the Cooperon 
mass in Eq. (13]): Dq 2 h-» Dq 2 + D/L 2 El with D/L% = 
2(-iEcos6 + A sm6 + A r]cos26). Thus, in the 2D ge- 
ometry, rj becomes a function of 9 and E, making the 
depairing term in the Usadel equation more complicated 
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than the simple AG term. However since the dependence 
of r\ on Le is logarithmically slow, we can evaluate it at 
the AG solution replacing 9 by 0ag(E). 



L, 



£1/2 



[24Ao(£' s - E)] 1 / 4 



E a 



E g -E 



(14) 



(we assume < 1). As a result, the depairing factor 
becomes energy-dependent: 



m _ K(T/T c ) ln min(L E ,L g ) 



9(9 ~ 9c) 



(15) 



where we had to introduce an infrared length scale L g to 
regularize the otherwise divergent rj(E — > E g ). Its ap- 
pearance is related to the breakdown of the mean-field 
approximation in the narrow region \E — E g \ < Ttaii 
[6j, where proliferation of instantons generates a finite 
correlation length L g ~ ^(r)(S 9 /rtaii) 1 ^ 4 - Substituting 
T) = v(E g ) into E^ G , we obtain for T = A(T) — E g : 



A(T) 



K(T/T C ) L g (T) 



In- 



45(3 - 5c) Co 
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(16) 



which in the zero-temperature limit reduces to Eq. ([T}. 

The theory developed above applies to quasi-2D films 
with the thickness d <C £o- For finite d < £0, there exists 
a contribution to r\ coming from short scales (I < r < d) 



where electron diffusion is 3D [241 ]: 



2A(T) 

TT 2 hD 
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Rq 



In 



rf 2 1 d 



(17) 



where p is the film resistivity, Rq = 27rfi,/e 2 , and £(T) is 
the temperature-dependent coherence length. Correction 
(JTTj which should be added to Eq. (fT5j) leads to a small 
increase of the width T. 

Subgap states. — A hard gap in the excitation spectrum 
predicted by the AG theory is smeared by disorder lead- 
ing to the formation of states at E < E g . These localized 
states are identified as instantons in the nonlinear equa- 
tion (UJ 0, (Hill [23. In Ref. 0, Lifshitz-type argu- 
ments were used to determine an optimal fluctuation of 
A(r) in Eq. §FZ§ giving rise to a finite DOS at E < E g . 
On the contrary, in Ref. [2l[, Eq. (fT2"j) was considered 
for a uniform A(r) = Ao, and instantons related with 
intrinsic nonlinearity of the Usadel equation were ana- 
lyzed. The results of Refs. Q and 21 1 are different and 



represent two asymptotics of a unique function of E g — E 
(a detailed discussion will be given elsewhere [26|). For 
small E g — E, nonlinearity is weak and the subgap DOS 
is due to optimal fluctuations of A(r) 0, whereas the 
mechanism of Ref. [2 11 ] is applicable only for very large 
Eg — E, where the DOS is exponentially small. 

Generalizing the 3D analysis of Ref. [6[ to the 2D case 
(24I . we find that the DOS decays exponentially in the 



subgap region: 



(p(E)) cx exp{ 

Itai] 

Mf) = 



[(E g - £)/Itail 

A(T/T c ) l2/3 



3/2 



0.47 



9(9 ~ 9c) 



(18) 
(19) 



[in the d-dimensional space, ]n(p(E)) cx — (E g — E) 2 ~ d ' 4 ]. 
At T = 0, Eq. (JTSJ reduces to Eq. Equation (JT5J) is 
valid as long as E g — E > r ta ii- Note that the instanton 
action obtained in Ref. [2l[ neglecting A(r) fluctuations 
is extremely large, S ~ g, already at E g — E ~ r ta ii- 
Hence, (p(E)) follows Eq. (fT8|) for all conceivable E < E g . 

Discussion. — Superconducting samples always have 
some amount of disorder which leads to smearing of the 
BCS density of states. We have considered the case 
of a minimal possible disorder — intrinsic randomness 
in a homogeneously disordered film responsible for the 
UCF in the normal state. The average DOS sketched in 
Fig. [1] is characterized by two energy scales. The shift 
of the gap edge Ao — > E g and related suppression of 
the CP height is controlled by the parameter T, so that 
Pma,x(E)/po w (A/r) 1 / 2 . At the same time, the width 
of the subgap tail is determined by a different parameter 
rtaii < r. Both r and rtaii are small in a clean system 
but get enhanced as the film becomes less conductive ap- 
proaching the critical point, g — > g c . Smearing of the 
DOS structure becomes very strong at g — g c ~ 1/.<? C , 
where our theory becomes inapplicable. The CP broad- 
ening is temperature dependent and grows at T — > T c 
due to the growth of the function K(T/T C ) (see Fig. [Z§. 

Our results for the disorder-induced DOS smearing 
should be compared with the smearing due to inelastic 
scattering In a 2D system, the inelastic rate due to 
Coulomb interaction is of the order of T ee ~ T/g at T ~ 
T c , and gets exponentially suppressed, T ee cx e _A ^ T ^ T , 
at T <C T c fioj ]. Therefore our mechanism always dom- 
inates at low temperatures. It also always dominates 
close to T c , where K(T/T C ) cx (1 - T/T c )- 3 / 2 . The 
general relation between the rates T and T ee depends 
on the proximity to the quantum critical point. Rela- 



tively far from it, at g 



> 



yfg~ c , inelastic scattering 



is the leading source of the DOS smearing at intermedi- 
ate temperatures, T\ < T < T2, where T± ~ Ao/ln/3 and 
T c -T 2 - P' 2 , with P = (g~g c ) 2/3 /g 1/3 - For films closer 
to criticality, g — g c \[g~ , disorder- induced smearing al- 
ways dominates over the inelastic smearing. 

Developed theory is expected to be most appropriate 
for very thin films of amorphous metallic superconductors 
(e.g., Mo-Ge, Nb-Si, W-Re, NbaGe), where suppression 
of T c by disorder is described by the fermionic mecha- 
nism [3, and R u > 1-2 kf2. 

Finally, we emphasize that the widely used phe- 
nomenological Dynes ansatz, p(E) = po Rc(E — iT)/[(E— 

ir) 2 -A 2 ]!/2 3 

is inapplicable when disorder is the main 
source of the broadening. The actual DOS profile then 
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depends on two energy parameters, T and r ta ii, and de- 
cays exponentially rather than algebraically at E < E g . 
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I. INVERSION OF THE FLUCTUATION PROPAGATOR 

Here we invert the matrix fluctuation propagator denned by the left-hand side of Eq. $7§: 

A^e) - KT^\{e,e')?^^k x {e') = <f>(e). (S20) 

Evaluating the derivative of F(e') = A(e')/[e' 2 + A 2 (e')] 1 / 2 , we get 

A l(£ ) - ^Y,M^ } F 0{ e')^l + [e t'"2g)]^ Al(£,) " m <S21) 

Since both A(e') and Ai(e') are logarithmically slow functions of e' , the second sum can be easily evaluated and we 
arrive at 

K 1 (e)-^TY J Ke,e')F (e')^^+L^\(e 1 T c )K 1 (T c )=4 > {e) 1 (S22) 
where Lq is the fluctuation propagator at zero frequency and momentum in the BCS theory: 

>2| 



The value of Ai(T c ) can be easily obtained from Eq. (|S22[) . In order to do this we multiply it by F(e) and sum 
over e. Using the SCE ([4]), wc immediately see that the first two terms in Eq. (|S22[) cancel and we obtain 

A 1 (T c ) = ^7rT^F o ( e )0( e ), (S24) 
where we use that, according to Eq. ([6]), A(T C ) = A(T). 

II. CORRELATION FUNCTION (AjAj) DUE TO MESOSCOPIC FLUCTUATIONS 

In this Section we evaluate the zero-momentum correlation function 

$ = (AimiAxCT^ (S25) 
A(T C )A(T C ) 

due to mesoscopic fluctuations of F^is and Xdis- 

Correlation function {Fdi s Fdi s ) 

The correlator of Gorkov functions in the Matsubara representation is calculated with the help of imaginary-time 
replica sigma-model following the line of Ref. [SI]. The resulting expression has the form 

(F dis (e, r)F dia (e', r')) = F (s)F (e') ^/^f )] ' , (S26) 

where v is the 2D one-particle DOS at the Fermi level (per single spin projection), and IT is the diffusion operator on 
top of the superconducting state: 

n"! =~DV 2 + €(e) + €(e'), (S27) 
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where 



For the zero Fourier component we get 

(F dis (e)F dis ( £ '))q=o = 



g(e) = yJe 2 + A 2 (e). 

F (e)F (e') 



4ir 3 v 2 D (£(e) + &(e') ' 
The corresponding contribution to Eq. (|S25[) has the form 

^(FF) Ll (2^T) 4 x ^ F ( £l )Fo (£ 2 )F (£3)^0 (64) 



ei,£2,e3,E4>0 



C(e a ) + C(e 4 ) 



(S28) 



(S29) 



(S30) 



This expression is typical to fluctuation contributions to $. Among four energy summations, two are logarithmic 
involving large energies, e 3> T c , while the other two come from e ~ T c , where A(e) can be approximated by A(T). 
The latter summations introduce the dimcnsionless function 



T, 



A(T) 



2>0 



£i€ 2 (€i + (£2) tt Jo cosh 



^ , 2 

■ tann 



A(T) 



2T 



cosh 6* 



14C(3) A(T) 
2tt 3 ; 



T c - T < T c 
T < T c . 



Performing summations over £2 and £4 in Eq. (|S30[) we get 
$(^) 



16tw 2 L>A 

Summation is done with the help of the SCE (|4]), and we obtain finally 



J' 2 



167rz/ 2 DA(r) u \T C 



T 



N 



(S31) 



(S32) 



(S33) 



Correlation function (5Adi s 5Adis) 

Mesoscopic fluctuations of the return probability which determines SX(e,e') = — A 2 ln[l/ max(e, s')t] have been 
calculated for the normal state in Rcf. [S2]: 

n c- u ^ r w 6\(e 1 ,e 2 )6\(e 3 ,e 4 ) ( 1 1 1 1 \ 

(A d i s (.£i,£2jAdi s ( l £3,e4j; q =o - ir , 2n 1 1 + 1 1 + I 1 + I 1 1 (oo4) 

where E»j = £, + £j . 

Generalization to the superconducting case is achieved by replacing e — ^ 2(e)- All four terms in Eq. (|S34[) equally 
contribute to and we get 

(AA) _ Ll (2nTr v /o( £ i)Fo(£ 2 )F (£ 3 )F (£ 4 ) 
* -Anr)4^^ ^ <5A(£ 1 ,£ 2 )a(£ 3 ,£ 4 ) g (£2) + g(e4) ' ( S35 ) 

Analogously to Eq. (|S30|) ]. summations over £2 and £4 yields N(T/T C ), while summations over £1 and £3 are converted 
to a logarithmic integral: 

^ - wJdMT) Ll (I) N (I) ^«coshA 9 ( C ,-0^ , (S36) 



and hence 

16tw 2 £>A(T) °V7 1 c y' V r < 



^ = 1fi ,./ nA ^ ^o ( J) * (£) [coshA 9 C* - I] 2 - (S37) 
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Correlation function (-Fdis^Adis) 
The cross term is evaluated analogously: 



^ FX) = ~2 ig 9 " ( N ( -J ) [coshA ff C* - 1]. (S38) 



16wv 2 DA(T) \T C J \T c 

Resulting expression for (AiAi) 
Adding (fS33|) . ([S37]) and (jS38|) . we obtain 

where we have introduced K(T/T C ) = Lq(T/T c )N(T/T c ), used the relation g = AttvD, and employed cosh 2 A g (* = 
.9/(3 — 5c)- With the help of Eq. (|S25j) . one readily obtains Eq. (fT0|) of the main text. 

It's worth noting that Eq. (|S39)) agrees with our previous results [S2]. Indeed, in the limit T — > T c , Eq. (|S39|) can 
be simplified with the help of Eqs. (|S23[) and (|S31j) as 

2T 3 

((A(T)A(T))) q=0 = K[i)v2D ^ {T) cosh 2 A sC *. (S40) 

On the other hand, the same correlator can be obtained from the correlation function of the coefficient a in the 
Ginzburg-Landau (GL) expansion [S2]: 

«aa» q=0 = cosh 2 A 9 C* (S41) 

with the help of the relation 

«A(T)A(T)» q=0 = (S42) 

where (3 = 7C(3)^/(87r 2 T 2 ) is the nonlinear coefficient in the GL functional. One can easily verify that Eq. (IS42[) 
coincides with Eq. (|S40[) . 

III. DENSITY OF THE SUBGAP STATES 

Introducing 8 = n/2 + itp, we rewrite Eq. (|12l) as 

t2v72„,, , _ <5A(r)sinh^ 



£'W + Fty) = -, (S43) 



where £ 2 = D/2A , and 



E 

F(ip) = — — — cosh 7/; + sinh?/> — f) sinh^ cosh-0. (S44) 
Ao 

At the minigap [E = E g ), cosh-0 g = yy^ 1 / 3 . For small deviation from the gap, the function F(ip) in the vicinity of 
its maximum can be written as 

Fty) w n(i> - fa) - 9 {t\> - V^o) 2 , (S45) 

where 



= ^^1-^/3./^—^, p=| ?? 1 /3 v / 1 _ ?7 2/3. (S46) 

A 2 
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Comparison of the linear in ip terms in the left-hand side of Eq. (|S43|) defines the relevant length scale 

L. - « . (S47) 



At the mean-field level, ip is real below the gap. Finite DOS corresponds to appearance of a nonzero Imip due to a 
large negative fluctuation of <5A(r). Its probability is given by 

V cx cxp (- J 5 A 2 (r) d d r^j = e~ s « , (S48) 

where we have used that the instanton size, Le, exceeds the correlation length, £oj of the order parameter fluctuations. 
At the quantitative level, the instanton action S can be estimated as follows. To produce a nonzero DOS at E < E g , 
the optimal fluctuation of <5A(r) should have the magnitude of — (E g — E) and the spacial extent of Le, which 
immediately gives the estimate [S3] 

Alt* (E g -E\ 2 - d ' A 

Sd W) \—) ■ (S49) 

To find the numerical coefficient in Eq. (|S49p . one has to solve the instanton equation. Measuring coordinates in 
terms of Le introduced in Eq. (|S47p . appropriately rescaling tp — ipo, and replacing sinhi/) in the right-hand side of 
Eq. (|S43|) by its value at E g , we rewrite Eq. (|S43[) in a universal dimensionless form: 

- V 2 4> + 4>- 4> 2 =h(r), (S50) 

where 

SA(r) = —^-h(r/L E ). (S51) 
p sinn ipg 

Minimization of the functional J h 2 (r) d d r leads to the fourth order differential equation for 4>(r): 

(-V 2 + 1 - 20)(-V 2 + (/)- 4> 2 ) = 0. (S52) 

The spherically symmetric optimal fluctuation solving (|S52[) in d dimensions satisfies the second-order differential 
equation [S4] 

-V 2 _ 2 + 0-0 2 = O, (S53) 

where V 2 _ 2 = d 2 /dr 2 — (d — 3)r~ 1 d/dr is the radial part of the Laplace operator in the (d — 2)-dimcnsional space. 
The instanton is characterized by the number 



i/2 



48tt/5, d = 3, 



J ^"-[u; t-i (S54) 

where 03 follows from the exact solution ^3(7") = (3/2) cosh _2 (r/2) [S3], while a 2 is obtained by a numerical solution 
of Eq. (IS531) . 

Returning to the dimensional variables, we get for the instanton action in the limit ^< 1: 

^m^r"- <*»> 



6 rf / 4 /(0) v E, 



J y 



In the 2D case, 



leading to Eqs. (TI5J and ([T5]). 



4a 2 L> A f Eg — E 
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IV. ROLE OF A FINITE FILM THICKNESS 



In films with a finite thickness d < £o> there exists a contribution to the depairing parameter 77 in Eq. (|13[) coming 
from large wave vectors (qd 3> 1), where diffusion is three-dimensional: 

2 ff(d) d 3 q 1 f 1 ' 1 

[here A = A(T) is the temperature-dependent BCS order parameter]. In this region, Coulomb effects are weak and 
all complications related with the energy dependence of A and A can be neglected. Equation ((5J is then replaced by 
a simpler expression written for an arbitrary 3D wave vector: 

5A(T,q) =L 9 (^)27rT^F dis (e,q), (S58) 

C £>0 

where L q (T/T c ) is the BCS fluctuation propagator at finite momentum: 

T^(TlT\ q-tV g(gW + 2A 2 (T) 

L q (T/T e ) = 2nT ^ (S59) 

In the limit q^Q 3> 1, one recovers the known inverse logarithmic decay of the fluctuation propagator |S3| : 

L q (T/T c ) « & » 1. (S60) 

With the help of Eq. (|S26p the correlation function /(q) = (<5A5A) q can be written as 

f(a)-L 2 ( T ) ( TA Y V / ^ n ^( k + q/ 2 ) n ^(k~q/2) 

where 1/3 = z//<f is the 3D DOS at the Fermi level. Integrating over k with the help of the Feynman's trick we get 

1 T 2( T \ (TA\ 2 V- 1 dx 

/(q) ~8^ AYJ [~^dJ £i ^ >Q g( ei )€(e 2 ) 7 Wl-^ + ( £ ( £l ) + £( £2 ))/D]i/2- ( } 

In the limit q£o 3> 1, only large £1,2 3> T c with €(e) ~ |e| are important. Replacing summations by integrations 
and using Eq. (|S60| we find 

/(q)= 8^ L 'fe) (^d) X c Rl - *)<7 2 + (ei + e 2 )/D}^ ~ 8^ (^d) " (S63) 

The obtained 1/g dependence of the correlation function f(q) leads in Eq. (|S64I) to the logarithmic contribution 
from the region I <C r <C d: 

leading to Eq. ([XT]). 

For even thicker films with d 3> £ 0j the 2D contribution (| 15[) is small and Eq. (|S64[) with d replaced by £0 gives the 
leading contribution to the depairing rate. 
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